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We analytically calculate weighted azimuthal asymmetries in semi-inclusive lepton-nucleon deep- 
inelastic scattering and Drell-Yan processes, using transverse-momentum-dependent partonic densi- 
ties obtained in a diquark spectator model. We compare the asymmetries with available preliminary 
experimental data, in particular for the Collins and the Sivers effect. We make predictions for other 
cases of interest in running and planned experiments. 
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H '. I. INTRODUCTION 

' Azimuthal asymmetries correspond to cross-section modulations depending on the azimutal angles involved in the 
^ process. Most of the time they are also spin asymmetries, in the sense that these azimuthal modulations appear 
with opposite sign when the spin of one of the participating particles is reversed. They are essential tools to study 
' ^ '' partonic transverse- momentum distributions (TMDs), defined as probabilities to find inside a hadron a parton with 
f-^ ^ longitudinal momentum fraction x and transverse momentum Pt with respect to the direction of the parent hadron 
I momentum [l| . TMDs can be used to construct a three-dimensional picture of partons inside hadrons in momentum 

space and may be related to the orbital angular momentum of partons ■ 
^ ■ In a recent paper, we published our results for all leading-twist TMDs in the context of a spectator model of 
the nucleon, using scalar and axial- vector diquarks, and further distinguishing between different isospin projections 
{ud and uu) [lo| . The free parameters were fixed by reproducing the parametrizations of Ref. and |12j for the 
unpolarized and helicity distributions at the lowest available scale, respectively, the latter being assumed as the model 
scale. Nonvanishing odd structures with respect to naive time- reversal transformations (for brevity, T-odd TMDs), 
' were generated by approximating the gauge-link operator with a single gluon-exchange interaction, representing the 
; rescattering of the struck quark with the spectator diquark at leading order in the strong coupling constant as- 
^-H ' Compared to our original publication, in this work we have modified the value we chose for at the model scale: 
\ instead of an ad-hoc nonperturbative value, we have computed it using renormalization-group equations at leading- 
■ order (LO). This change has the effect of resizing all T-odd functions by a global constant. 

] A first test would be to calculate the Pt dependence of unpolarized cross sections, either in semi-inclusive deep- 
T— I . inelastic scattering (SIDIS) or in Drell-Yan hadronic collisions. However, asymmetries are usually preferable, because, 
^ ■ being defined as ratios of cross-section combinations, are less sensitive to systematic errors and to theoretical uncer- 
tainties as well. 

Two types of azimutal asymmetries are commonly considered: weighted and unweighted. From the experimental 
point of view, unweighted asymmetries are easier and safer to measure. There is by now a good amount of data, 
mainly on single-spin asymmetries (SSA) in SIDIS [Tsl. [l^. Experimental data for weighted SSA are scarce, with low 
statistics, and still preliminary .ISJ. However, from the theoretical side unweighted asymmetries are more complex 
because transverse momenta are intertwined in a convolution that can be broken only by assuming a specific 
distribution, typically Gaussian. Weighted asymmetries are preferable: they can be written in a model-independent 
way in terms of "collinear" objects — parton distribution functions (PDFs) or pj- moments of TMDs, and analogously 
for fragmentation functions Most of the time, this simplifies the theoretical treatment, allows a simpler study of 
factorization and scale evolution, and allows more freedom beyond the choice of a Gaussian ansatz for TMDs. 

In fact, our model TMDs display a Pt distribution which is not Gaussian (while also displaying a strong flavor 
dependence and strong x — p^ correlations) [1ft] : as such, it does not allow to analytically work out unweighted SSA as 
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simple products of terms, and requires a numerical approach [l7j (see also the approximations discussed in Ref. fisj). 
On the contrary, weighted asymmetries in our model can be always calculated analytically. 

For what concerns the scale evolution, equations are known only for the first Pt moment of the Sivers function 
function jT9l - l2]| , and have a non-diagonal form that makes their treament more difficult than for the standard collinear 
PDFs. In order to include some evolution effects while avoiding these complications, we compute the moment of the 
Sivers function at different scales using only the diagonal part of its evolution equations, and we extend this treatment 
also to the other TMD moments. 

For all the above reasons, we choose to calculate weighted azimuthal (spin) asymmetries and to compare with the 
few available data: the double-spin asymmetry All in SIDIS with longitudinally polarized protons [24I, the Collins 
and the Sivers effect in SIDIS 15], all measured by the HERMES collaboration. There are also All data collected 
by the HERMES and COMPASS collaborations with a deuteron target (see for example Ref. [1^), but in this paper 
we will consider only measurements directly on a proton target. We will show predictions, instead, for other cases of 
interest in view of running or future experiments. For example, at JLab the E06-010 experiment (23 | is measuring 
the Collins and the Sivers effects using the 6 GeV energy beam hitting on a transversely polarized ^He (effective 
neutron) target, in order to extract the neutron transversity distribution. As for Drell-Yan, the fully polarized 
process with antiprotons is planned to be measured at FAIR (GSI), in order to perform a self-consistent extraction of 
the transversity distribution [25l - [29| : while the COMPASS collaboration is planning to measure a SSA using a high 
energetic pion beam on a transversely polarized proton target at CERN [s^'ssl, where the T-odd Sivers function 
could be extracted and its predicted non- universal behaviour directly tested 34] . 

The paper is organized as follows. In Sec. [Hi the kinematics and the formulae for the weighted azimuthal (spin) 
asymmetries in SIDIS and Drell-Yan are briefly illustrated. In Sec. Illli the formalism and the analytic results for 
T-even and T-odd TMDs in the spectator diquark model of Ref. are recalled. In Sec. IIVI the updated model 
T-odd TMDs are compared with the available parametrizations and the weighted SSA are compared with the few 
SIDIS experimental data, making predictions for asymmetries in the Drell-Yan process and for other cases of interest. 



II. WEIGHTED AZIMUTHAL (SPIN) ASYMMETRIES 

In the following we will make use of light-cone coordinates. We introduce the light-like vectors n± satisfying 
nj_ = 0, n+ • n_ = 1, and we describe a generic 4- vector a — [a~ , a'^ , a^], where — a ■ Uzp. We will also make use 
of the transverse tensor — ^^'^''^n^^n^^, whose only nonvanishing components are e]? — — e^^ — 1. 

We consider the SIDIS process of a lepton on a (polarized) nucleon, as well as hadronic collisions leading to Drell-Yan 
lepton pairs. 



SIDIS 



In a SIDIS process a (polarized) lepton with momentum I is scattered to a final momentum V by a nucleon target 
with mass M, momentum P, and spin S, leaving also a final hadron semi-inclusively produced with mass Mh and 
momentum P^- The 4-momentum transferred is denoted hy q — I — I' , with = —q^. We introduce the invariants 
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where k = p + q, and p is the momentum of the parton before the scattering. We choose a set of frames where P has 
no transverse components, i.e.. 



P 



M2 



2P+ 



,P+,0 



(2) 



Hence, the target polarization 4-vector and the parton momentum can be parametrized as 

XhM XhP^ 



2P^ 



M 



P 



2xP+ 



, xP^ , Pt 



(3) 



with Xh the hadron helicity. Analogously, Ae represents the incoming lepton helicity. 
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In single-photon-exchange approximation, the cross section can be parametrized in terms of 18 structure func- 
tions [33 • Here we are concerned with the regime ^ P^_|_ , . An analysis based on TMD factorization reveals 
that only the following 8 structure functions are present at leading order in a M^/Q^ expansion (twist expansion) ^: 
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+ B{x,y) sm{3ct>,-cl>s)F^f'^'^-^'^ 



+ \S^\ Xe C{x, y) cos(0„ - cbs) F™^^'^'--^^) 



where a is the fine structure constant, s — Q^/xy is the center-of-mass energy squared, 7 — 2Mx/Q, 'pin'Psj 
the azimuthal orientations of Ph and S with respect to the scattering plane, respectively (see Fig. [TJ for their formal 
definition, see Eqs.(2.3)-(2.6) in Ref. [H), and 
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Note that the above expressions include mass corrections, related to 7, and the factor l/(a:^y^), which must not be 
dropped when integrating separately numerator and denominator of the asymmetries (see Sec. II V[) . 




In Eq. Q, each contribution F^y depends on x,z,P^j_, and on Q^; it refers to a lepton with polarization state 
X (that can be unpolarized, U, or longitudinally polarized, L, with helicity Ag), and to a target with polarization 
state Y (that can be unpolarized, U, longitudinally polarized, L, with helicity A/i, and transversely polarized, T, with 
polarization The superscript W refers to the azimuthal distribution of detected hadrons. At leading order in 



^ Since we are keeping only the leading-twist contribution, we have shortened the original notation Fm/ x, F^^l^'^ 't's ) ^ q£ Kef. [35(1 t° 
^UU : F'^^'^'^ 't's) ^ since we are not sensitive to longitudinally polarized virtual photons. 
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the strong coupling (LO) and leading twist, the structure functions read [3£ 
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where h = Ph^/\Ph±], and the convolution is defined by 
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with the sum running over all parton and antiparton flavors a. 

The different structure functions can be extracted thanks to their specific dependence on the azimuthal angles. 
Upon integration in dP^^, only the convolutions in Fjjif and Fll can be solved analytically, because of the simple 
weight w = 1. In this case, the well known double-spin asymmetry can be formed, 
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giving access to the helicity distribution gi{x) = J dp^ (7°^(x,pJ). This asymmetry is analyzed in Sec. IIVB] 

In the other cases, an analytical solution can be achieved only with a suitable model dependence upon p^ and 
fey of the distribution and fragmentation functions, respectively [l^. The most widely used Ansatz is a Gaussian 
distribution in p^ and k^ (usually also assumed to be flavor independent and independent of x and z). It allows 
an analytic calculation of all convolutions in Eq. (|9]). However, theoretical calculations indicate that perturbative 
corrections generate non-Gaussian tails of the TMDs when P'^^_ ^ M"^ [s^ and available experimental data poorly 
constrain the p^, dependence of partonic densities in general, leaving room for Ansatze different from the Gaussian one. 
Actually, our model displays a p^ distribution which is not factorized, has a strong flavor dependence, and above all 
is not Gaussian As such, it does not allow to analytically calculate all the convolutions in Eq. ([9]). In this paper 
we will consider only properly weighted asymmetries, which break the convolutions in a model-independent way 

and result in factorized expressions such as the one in Eq. (1101) . We leave the computation of unweighted asymmetries 
to a future work, where we will consider a direct (numerical) calculation of all the convolutions needed in Eq. ([5]) [l3l ■ 
Weighted asymmetries are defined as 
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where daxY refers to the contribution in Eq. Q for the lepton probe with polarization X and nucleon target with 
polarization Y . Typically, the weight W can be function of 4>s, 4>h, and of suitable powers of Qt = \Phi.\l z. 
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For SIDIS, the following weighted azimuthal asymmetries are relevant and theoretically clean [37[: 
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/(")(:.) = J dp^^^j fix,pl) (15) 

is the n-th pT-nioment of a parton distribution / (and analogously for a fragmentation function). 

In every parton density above (or p^ moment of it), the dependence on the hard scale is understood. Evolution 
effects with running are well known for /i, gi, hi, and Di. From the discussion in Sec. 7.3 of Ref. [SG*] and the 
explicit study of Ref. [20|, it should be evident why we can hope to study the weighted asymmetries (I12l - I14p in a way 
similar to the double-spin asymmetry of Eq. ([TU)) . including the scale dependence of the involved functions. Explicit 

calculations are available only for the evolution equations of the Sivers function f^r^^^ [H-ill- They are nondiagonal 
and thus considerably more complex than those of twist-2 coUinear PDFs. However, their diagonal part is identical 
to the evolution equations for /i and should be dominant at high x [20| . Calculations of the real perturbative 
contributions in giT [ssj suggest that the evolution equations for this function will have a diagonal part identical to 
that of gi. In summary, we approximately implement the effect of scale evolution in the weighted asymmetries (jl2l 
- UM by evolving f^r^^\x) and g[^{x) using the same LO evolution of fi{x) and gi{x) (see, e.g., Ref. 39]). For 

H^'"^\z), we used the same LO evolution as the transversity fragmentation function Hi{z) (see, e.g., Ref. 40]). The 
asymmetries (fT^ - [Hl) are considered in Sec. IIVB] 

In our analysis, we will not take into consideration the following asymmetries 
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because the higher powers of in the weight emphasize the high-p,, components of the asymmetry, that are 
dominated by perturbative QCD corrections [3a] and do not give enough information about TMDs. In the above 
equations, we nevertheless gave the explicit expressions for these weighted asymmetries in the TMD framework, 
since it might still be possible to isolate their TMD component by means of differences, or ratios between different 
measurements, that cancel the perturbative contributions. We will not discuss these possibilities in the present work. 



B. Drell-Yan 



In a Drell-Yan process, a pair with individual momenta ki and fc2, formed by a lepton and its antilepton partner, is 
produced from the collision of two hadrons with momentum P^, mass Mj, and spin Si, with i = 1,2. The cm square 
energy available is s = (Pi + P2)^; the momentum transfer is now time-like and gives the invariant mass of the pair, 
i.e. 



g2 = M2 = (fci-ffc2)2 = (pi+p2)2, 



(20) 
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where pi are the momenta of the annihilating partons, with i = 1,2. If and are the dominant hght-cone 
components of hadron momenta in the regime where the factorization theorem holds |41| , we can define the following 
invariants: 
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We choose Xi as the momentum fraction of the parton in the beam; namely, 
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In the following, we will plot SSA as functions of a:^? and/or y at a given r: this choice will probe very different regions 
in xi and X2, where the model parton distributions involved behave very differently. And this feature will show up 
clearly in the asymmetry. 

Because of momentum conservation, the 4-momentum transfer can be parametrized as 



X2P2 ,XiP^,PiT +P2J 



(23) 



If the transverse momentum of the final lepton pair is Qt = Pit + P2t 7^ 0, the directions of the two annihilating 
partons are not known. Hence, it is convenient to select the socalled CoUins-Soper frame [i^ described in Fig. [2] 
The final lepton pair is detected in the solid angle {9,(p), where, in particular, (j) (and all other azimuthal angles) is 
measured in a plane perpendicular to the indicated lepton plane but containing h = Qt/IQtI = Qr/qT- 




FIG. 2: The CoUins-Soper frame. 



The fully differential LO cross section can contain many contributions according to the polarization state of the 
two colliding hadrons ^43]. Having in mind the phenomenology that could be explored at GSI-FAIR and CERN, 
or at J-PARC and RHIC, we limit the physics case to the most relevant combinations. The master channel is the 
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rl X reaction [25l - l27f . whose polarized g^— integrated cross section is simply [44 
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where dfJ = sin 9d9d4> and IS'^tI, are the transverse polarization and the azimuthal angle of its vector for hadron 
i. Combining Eq. (|24p with the unpolarized part, we can form the well known double-spin asymmetry 
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which gives direct access to (the valence part of) the transversity distribution hi . Monte Carlo simulations about the 
feasibility of this measurement were presented in Ref. [2^ . Predictions based on our diquark spectator model will be 
shown in Sec. IIVCI 

When in the collision only one hadron can be transversely polarized, the fully differential leading-twist cross section 
reads [H,!!! 
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The structure functions F^y depend on xi,a;2,qT- They read [i^ 
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The convolution C is defined as 
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As in the SIDIS case, only the convolution in Fjju can be solved analytically upon integration in dq^, because of 
the simple weight w = 1. Again, instead of introducing a suitable model dependence upon transverse momenta, we 
follow the strategy of considering properly weighted asymmetries, now defined as 
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where daxY refers to the contribution in Eq. (j26p for the hadrons 1 and 2 with polarizations X and Y, respectively. 
Typically, the weight W can be function of (t>s2i 0: ^-iid of suitable powers of q^- 
We consider the following weighted asymmetries: 
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For the same considerations expressed in the case of SIDIS, we refrain ourselves from considering the following 
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asymmetries 



where the right-hand-side represents only the TMD contribution to the asymmetry, which is suppressed with respect 
to perturbative contributions. 

Monte Carlo simulations about the feasibility of measurements of the asymmetries (PT|) and ([5^ at GSTFAIR, 
COMPASS, and RHIC, were presented in Refs. [3l|, IH, H^, • Predictions based on our diquark spectator model are 
displayed in Sec. IIV CI Similarly, while in Eqs. pil32l) the dependence on the hard scale is understood, evolution 
effects with running are included in the same way as discussed at the end of the previous Sec. Ill Al The chiral-odd 
function /i^j^^^^ is assumed to follow the same LO evolution equations as the chiral-odd function hi. 



III. T-EVEN AND T-ODD PARTONIC DENSITIES IN A DIQUARK SPECTATOR MODEL 

The TMD showing up in Eqs. (|12m4l) and pil32p can be calculated analytically in the spectator diquark model. 
They can be extracted from the basic quantity named quark-quark correlator [sBj , i.e. 

^{x,p^,S)^ I e^P< {P,S\mi^Mm\P,S)\^^^^, (35) 

where 

U[o^^] = V e-'3 ^0 d^-M^) (36) 

is the socalled gauge link operator connecting the two different space-time points and ^ by a specific path followed 
by the gluon field A, which couples to the quark field ■0 through the coupling g. The gauge link ensures that the 
matrix element of Eq. psp is color-gauge invariant. 



P 
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FIG. 3: Tree-level cut diagram for the calculation of T-even leading-twist parton densities. The dashed line indicates both 
scalar and axial- vector diquarks. 

The correlator ([55]) can be evaluated in the spectator approximation [l^ : a completeness relation is inserted 
and at tree-level the sum over final states is truncated to a single on-shell spectator state with mass Mx (see Fig. [3]) , 

cl>(x,p.,5)^-l^ ^ M^'\S)M^'\S) (37) 
(27r)-^ 2(1 - a:)P+ p^=t(x,p^) 

where the on-shell condition {P — pY — Mj^ for the spectator implies for the quark the off-shell condition 

p^ = T(x,p^) = -^-^-^^^^ + m^, Lj,im'')^xAa + il-x)m'' -x{l-x)M^ . (38) 

1 — a; 

We assume the spectator to be point-like, with the quantum numbers of a diquark. Hence, the proton can couple 
to a quark and to a spectator diquark with spin (scalar X = s) or spin 1 (axial- vector X = a), as well as with 
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isospin projection (wd-like system) or 1 (ww-like system). Therefore, the tree- level "scattering amplitude" Ai^^^ is 
given by (see Fig. [3]) [13] 



2 

■ ys U{P, S) scalar diquark, 



X(°)(5) = (P-p|^(0)|P,5)= <; ^^7"^ (39) 

— ; e*(P — p, Aq) U (P, 5*) axial- vector diquark, 

f — m 

where the nucleon-quark-diquark vertex can have the form 

a iv^) 

ys = igsip^) 1 scalar , = i '^^y' 7^ 75 axial-vector , (40) 

v2 

with gx (p^) a suitable form factor. The e^(P— p, Ag) is the 4- vector polarization of the spin-1 axial-vector diquark with 
momentum P — p and helicity states Aq. In Ref. [lO|, several choices were analyzed for the diquark polarization sum 
^fiiy _ ^ ^ and for the form factor gx{p'^), with all analytic formulae listed in the appendices. Numerical 

results were presented only for d'^'^ propagating transverse diquark polarizations, and for the so-called "dipolar" choice 

gx{p)^9x— -^^9x 2 ■ (41) 

Here, we keep the same choices and in the next subsections we list the results for the T-even and T-odd leading-twist 
TMDs. 

A. T-even 

The prototype of TMDs is the unpolarized parton distribution, which can be defined in the following way: 
/i(x,p,) = i Tr [(<i>(x,p,, S) + $(x,p,, -5)) 7+] + h.c. 

= i 2(1 - .)P+ [{m''\s)M^'\S)+m'''\^S)M^^\-S)) 7^^ ^"'^ 



h.c. 



and can be computed by inserting in A^'"' of Eq. (15^)) the rules (PH)) for the nucleon-quark-diquark vertex, the dipolar 
form factor of Eq. (1411) , and the diquark polarization sum. 

To express in a condensed way all the other TMDs, it is convenient to introduce the two combinations 



/i (x,pj_2(^/i (x,p,) gii(2;,p,)J- 2f«2+L2^A2)14 ■ (4b) 



5? 


(m -1- a;Af)2 (1 - xf 


(27r)3 


2[p2+L2(A2)]4 


.9? 


p2 (1 _ xf 


(27r)3 


2[p2+i2(A2)]4 ' 


5a 


p2 (H-a;2)(l-a;) 


(27r)3 


2[p2+L2(A2)]4 ' 




(m -f a;M)2 (1 - a;)3 


(27r)3 


2[p2+L2(A2)]4 



Then, the TMDs can be written as (compare with Ref. flo| ^ 



!r\x,Pl) = fr'\x,pl)+f,-'^^'>{x,pl) , (47) 

ff''\x,pl) = It'^''\x,pl) + f;'''-''\x,pl) , (48) 

gfl\x^pl) = ft'^'\x,pl) - fr^'\x,pl) , (49) 

gtr\x,pl) = ft'^^\x,pl) ~ fr^''\x,pl) , (50) 
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q(s 
9lT 



q(a 
9lT 



IL 

± q{a 
IL 



'IT 
'IT 



'it 

±q{a 

'it 



hi 



q{s 



x,pI) = 



2M 



m - 



— X m + xM 



x,pI) = , 



2 N 



5iT (^^'Pt) > 



x,pI) 
x,pI) 



ff'Hx,pl) 



2x 



2M2 



hp{x,pl). 



■f^'^'\x,pl), 



(51) 
(52) 

(53) 
(54) 

(55) 
(56) 

(57) 
(58) 

(59) 
(60) 



We stress that the previous relations among TMDs should be considered simply as a convenient way to write the 
results of our model, and should not be considered as general, as can be also deduced by the fact that they are all 
different in the scalar and axial- vector diquark case. For instance, the relation involving the pretzelosity h^j,^ recently 
discussed by several authors J,. i49i-i5l]. holds in the scalar diquark case, i.e. 



dld^^PT) - hl{x,pj 



P. 



2M2 



hiTix,PT) ; 



(61) 



but it is not true for the axial- vector diquark, as it can be deduced from Eqs. ([5(I)l . ([SO)) , and ([55)1 . by noting that 
h^rf- — (see also Ref. fiol)- Therefore, such a relation cannot be totally general. 

The Pr-integrated results lO] can be expressed in terms of the pr-integrated combinations (1431146^ : 



/+'(^)(x) 



gl 2{'m + xMf [1 - xf 



^2 



9 s 



24L6(A2) 

{l-xf 



(2^)2 24Lf(A2) 



(62) 
(63) 



The final result is 



ft''^^\x) 
f-'^^\x) 



gl il + x^){l-x) 



(2^)2 24L4(A2) ' 
gl 2{m + xMf (1 - xf 



{2ny 



24L6(A2) 



ff'\x)^f+''^^\x) + f-''^^\x), 
ft''\x)^ft'^''\x) + f;'^'''\x), 



(64) 

(65) 
(66) 

(67) 
(68) 



gf\x) = ft'^'\x)-f;''^'\x), 
9t\x) = ft'^''\x)^f^''^''\x)^ 



(69) 
(70) 
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2x 



1+^2 

The weighted SSA of Eqs. (IT4| . pT)) . (fTS)) . and (p4)) . require the knowledge also of the foUowing p^-nionients: 

gl (m + xAf) (1 - xf 



(2^)2 24ML4(A2) ' 
5^ (m + xM) X (1 - a:)2 
(2^ 24ML4(A2) 



(71) 
(72) 

(73) 
(74) 



/j^/(«)W(a;) = ^5^^"^^^^(x) 



9(a)(1), 



(75) 
(76) 



[l-xf 



(27r)2 24M2L2(A2) 



(77) 
(78) 



For completeness, we also list the py-moments of the combinations (j43ll46p . which are useful to calculate the mean 
squared transverse momentum associated to each of the above TMDs: 



^+.(.)(i)(^) 



^2 (m + a;Af )2 (1 - xf 



(27r)2 96M2L4(A2) 

gl (m + a:M)2 (1 - xf 
(2^ 96M2L4(A2) 



(79) 
(80) 



/i 



-g(a)(l) 



(^) 



5? 2 (1 - xf 



(27r)2 96M2L2(A2) ' 

gl 2{l + x^){l-x) 
(27r)2 96M2L2(A2) 



(81) 
(82) 



p- P 



FIG. 4: Interference between the one-gluon exchange diagram in eikonal approximation and the tree level diagram in the 
spectator model. The Hermitean conjugate diagram is not shown. 



B. T-odd 



There are only two T-odd structures that can be extracted at leading twist from the quark-quark correlator of 
Eq. ([55)1 . namely 



'-^Ip^ fM^,pl) = ~l Tr [mx,p^,S) - <i>(x,p., -5)) 7+] + h.c. 



(83) 
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*J 1 

hi{x,pl) = - Tr 5) + <i>(x,p,, -S)) ia'+ 75] + h.c. , (84) 

which are usually referred to as the Sivers [H^ and Boer-Mulders [13] functions, respectively. 

Using the diagram rules for A^^") in Fig. [31 from Eq. ([57)) we find that these expressions vanish at tree level, 
because there is no interference between two competing channels producing the complex amplitude whose imaginary 
part gives the T-odd contribution. We can generate such structures by considering the interference between the tree- 
level scattering amplitude A^'"' and the single-gluon-exchange scattering amplitude A^*-^^ in eikonal approximation, 
as shown in Fig. [31 (the Hermitean conjugate partner must also be considered). This corresponds to the one-gluon- 
exchange approximation of the gauge link operator of Eq. p6p [SSJ . 

The Sivers and Boer-Mulders functions can then be rewritten as 

%^Ai(.,p^) = -i^^^^-i^Tr [[m^'\S)m''\s)-M^'\-S)m''\^S)) 7+] + h.c. , (85) 

"-^hiix^pD - ^ 2(1 - x)P+ [[m^'\S)M^"\s)+M^'\~S)M^'\-S)) ^a»+75] + h.c. . (86) 

Following the lines of Ref. flo| . we come to the final form with scalar and axial vector diquarks: 



ht''^'\x,pl) = f^^'^^^\x,pl), (89) 
/i^'(")(x,p2) = _i/^^/'^)(^,p2)^ (go) 

with Cf = 4/3. 

The high— Pt tail of above expressions could break the positivity bounds [H^l- However, the spectator model is 
supposed to be valid for p^. ~ and for reasonable choices of the parameters no problems with positivity occurr in 
this region (for more details, see Ref. [l0|V 

Finally, in the weighted SSA we need the first moments of such functions, which read [Toj 

^±q(s){i}^ ^_ 9s TrCpas im + xM){l~xf 
8 (27r)3M K(A2)]2 ' 

r±q{a){l)f ^ _ 9I TT CfQ^ X {ui + xM) {1 - x)^ 

^^^-8(27r)3M [L2(A2)]2 ' ^^^^ 

^^^-L^9(«)(l)(^) , (93) 

/i^«(«)(i)(x) = -i/^-L/-)(i)(^) . (94) 

IV. RESULTS AND PREDICTIONS 

In the following, we present the results for the weighted SSA discussed in Sec. [ill We recall that the model 
parameters were fixed in Ref. [l3| by fitting the ZEUS parametrization for and ff at the lowest available scale 
Ql = 0.3 GeV^ [11], and the GRSV2000 parametrization at LO for gf and gf at = 0.26 GeV^ [Tl|. The connection 
with the model TMDs is realized by 

fu _ 2 ru{s) 2 ru{a) 
Jl — 7l + J-^ 

ft = c^ft^'''\ (95) 

with the coefficients c^' given in Tab. I of Ref. [l^. They give the relative probability weight of each quark-diquark 
configuration in spin and isospin space: scalar-isoscalar diquark s with quark u (c^), axial vector-isovector diquark 
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o with isospin projection and quark u (ca), axial vector-isovector diquark a' with isospin projection 1 and quark d 
{<)■ 

The model (hadronic) scale is assumed to be Qq = 0.3 GeV^. At this scale, the contribution of sea over valence 
quarks is approximately within 10-20% for x > 0.1; therefore, we have involved only valence contributions in the 
fitting procedure, keeping in mind that there will be discrepancies at x < 0.1. Viceversa, for the Di fragmentation 
function we employ the parametrization of Ref. [ssl j where both favoured and unfavoured channels are included, the 
latter being particularly important even at low Q^. For the Collins function in Eq. (fT2l). we adopt the analytic 
expression of Ref. ^6], which was obtained in a similar spectator approach as our model TMDs. However, in this case 
the unfavoured contributions are simply assumed to be equal and opposite to the favoured ones. 

Evolution with running is implemented by suitably extending the HOPPET code of Ref. [s^] to include also 
chiral-odd partonic densities. The modification is here considered at LO only, and the evolution of all partonic densities 
will be considered at LO as well in the following, because the factorized expressions of the SSA in Eqs. (|12p - p^ 
and (I3T]),(I32]) are vahd at LO only. 

In the following subsection, we reconsider evolution properties of our T-odd TMDs. In the next two subsections, 
we will present our results for weighted SSA in the SIDIS and in the (polarized) Drell-Yan processes, respectively. 

A. Results for model T-odd TMDs 

As already anticipated in Sec. Ill A| evolution equations are not available for all partonic densities in the equations 
for the considered SSA. At LO, the DGLAP equations for f^rp^^\x), g[l^{x), hf'"^\x), are assumed to be the same as 
for /i(x), gi{x), hi{x), respectively. 

Moreover, the first moments of these T-odd TMDs depend linearly on Us [see Eqs. ((9T |) -(|94 |) ]. Therefore, a 
consistent treatment of evolution effects requires to deduce from the LO renormalization group equations the value 
of as at the scale of the hadronic model, Qq, while asiQ^) is often considered as a free parameter. 

In this sense, in the following we update the results for the model T-odd TMDs with respect to our previous 
paper [l^l by replacing the ad-hoc parameter asiQo) = 0-3 with the value of as{Qo) — 0.697 deduced from the LO 
running of the strong coupling starting from a value of a(M|) = 0.125 (see discussion in Ref. and using the 
standard choices of the HOPPET evolution program (57| . 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

X X 



FIG. 5: The first Pt moment x/^y [x) of the Sivers function; left (right) panel for up (down) quark. Darker shaded area: 
statistical uncertainty band of the parametrizations from Ref. l59l |. lighter shaded area from Ref. [60]. Solid line for the results 
of the spectator diquark model at the parametrization scale Q = 2.5 GeV^. 

In Fig. [51 the first moment xf^r^^^ {x) of the Sivers function is displayed for the up and down valence quarks in the 
left and right panels, respectively. The darker shaded area represents the uncertainty due to the statistical error in the 
parametrization of Ref. 59] , while the lighter shaded area corresponds to the same in Ref. ^60>] . Both parametrizations 
are deduced by fitting data for the Sivers effect measured by the HERMES [6l| and COMPASS [6l,[63] collaborations, 
hence at an average scale = 2.5 GeV^. The solid line is given by combining the model results of Eqs. (|9T|) and ([92]) 
according to Eq. and by further evolving them at LO from the model scale Qq = 0.3 GeV^ to the above mentioned 
parametrization scale. 

First of all, we observe the agreement between the signs of the various flavor components, which also agree with 
the findings from calculations on the lattice 64^. The agreement between model and parametrizations is very good 
for the up quark, even if the maximum is slightly shifted towards higher x. The size of the function is instead too 
small for the down quark, and its shape shifted such that the maximum occurs at a higher value of x « 0.3. 
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FIG. 6: The modulus of the first p-r moment \xh'^'^^\x)\ of the Boer-Mulders function; left (right) panel for up (down) quark. 
Darker shaded area: statistical uncertainty band of the parametrization from Ref. [gSj, extracted from SIDIS data at = 2.5 
GeV^. Lighter shaded area: error assigned to parametrization from Ref. [6^ . extracted from Drell-Yan data at = 1 GeV^. 
Solid and dashed lines for the results of the spectator diquark model evolved at — 2.5 and 1 GeV^, respectively. We remark 
that the results of the spectator diquark model are negative for both up and down quarks (see Ref. [l3l); here, we plot their 
modulus in order to conveniently compare with the published parametrizations. 



In Fig. [SI the modulus of the first moment \xh-^ (2^)1 of tli6 Boer-Mulders function is displayed for the up and 
down valence quarks in the left and right panels, respectively. The darker shaded area represents the uncertainty due 
to the statistical error in the parametrization of Ref. [65j , deduced by fitting the cos 20 asymmetry in the unpolarized 
SIDIS data of Ref. at = 2.5 GeV^ [see Eq. dH])]. For the down quark, the uncertainty in the fitting 

parameter is so small that the area reduces to a thin line. The lighter shaded area corresponds to the similar analysis 
from Ref. ^6^] but for the unpolarized Drell-Yan data of Ref. 70] at = 1 GeV^ [see Eq. (|33p ]. The solid line is 
given by combining now the model results of Eqs. and (|M)) according to Eq. and by further evolving them 
at LO from the model scale Qq = 0.3 GeV^ to = 2.5 GeV^. The dashed line displays what we obtain when we 
perform the LO evolution up to = 1 GeV^. Hence, the solid line must be compared with the darker shaded area, 
while the dashed line with the lighter shaded area. 

We remark that the cos 20 asymmetry in Drell-Yan involves two Boer-Mulders functions and the fitting procedure 
cannot fix the sign of the function depending on the fiavor involved. The cos 20 asymmetry in SIDIS involves 
in combination with the unknown Collins function, whose sign in turn depends on the sign of the transversity hi. 
The extraction of Ref. j65j assumes the favoured Collins function to be positive and obtains negative Boer-Mulders 
functions for both u and d quarks, in agreement with our model TO], with lattice calculations |;64], and also with 
other models [tiI, [t^]. The size of \xh'^ ^^'^^ {x)\ is too high in Fig. [SI while the result for the down quark seems in 
better agreement with the SIDIS parametrization than with the Drell-Yan one. However, the comparison should be 
considered with some care because our model results contain only the pure contribution from valence quarks and the 
data are contaminated by perturbative QCD contributions and higher twists (e.g., the Cahn effect in SIDIS). 

B. The SSA in SIDIS 

In the following, all displayed experimental data for weighted SSA in SIDIS were collected by the HERMES 
collaboration at the cm energy squared s = 56.2 GeV^ and the average scale (Q^) = 2.5 GeV^ (but evolution effects 
were included for each Q'^{x) in the x distribution). Applied experimental cuts are: 

0.023 < a; < 0.4 , 0.2 < z < 0.8 , 

Vminix) < y < 0.85 , ymin{x) = max 

Particular care must be taken when considering the asymmetry as a function of z, namely when separately integrating 
its numerator and denominator as functions of x and y. 

In Fig. [3 we show the x dependence of the virtual photon asymmetry Ai from Eq. (jlOp . Experimental data for 
both 7r+ (left) and 7r~ emission (right) are taken from Ref. [l^ in the above mentioned SIDIS kinematics of HERMES. 

In all cases, the solid line is the result of our spectator model when employing the fragmentation function Di of 
Ref. [sll, including LO evolution of all partonic densities to the experimental scale. The agreement with data is 



/r2 



1.1, 



1 4 - M 

•:{s~NPy [l-x) {s- M^)_ 



(96) 
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FIG. 7: The double-spin asymmetries Ai{x) from Eq ^lOp (see text) in the SIDIS kinematics of HERMES with emission of tt^ 
(left) and -k~ (right). Experimental data from Ref. [2^. Solid line represents the result of the spectator model. 



satisfactory. The deviations from data at low x are driven by contributions of nonvalence partons (sea quarks) , which 
are not included in the present version of the model (see the comment at the beginning of Sec. lIVp. 
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FIG. 8: The weighted single-spin asymmetry A'^'^ ^M'kh+'I's) in the SIDIS kinematics of HERMES (Collins effect) 

with emission of tt^ (left) and vr" (right), as a function of x (above) and z (below). Experimental data from Ref. ^15,]. Dashed 
line for the result of the spectator model at its scale Qq = 0.3 GeV^, solid line for the result at the experimental scale — 2.5 
GeV^ (see text for details about evolution). 

In Fig. [SI wc show the x (top panels) and z (bottom panels) dependences of the weighted SSA A'^'^ sin{(t,h+(ps) £j.qj-^ 
Eq. for both tt"*" (left) and tt~ emission (right) in the same SIDIS kinematics of HERMES as in the previous 
figure (Collins effect [231 )■ Experimental data are taken from Ref. [15J- The dashed line represents the result of the 
SSA when calculating it at the model scale Qq = 0.3 GeV^, where the analytic expression of the Collins function 
is taken from the consistent spectator approach of Ref. f56j; the Di of Ref. [53] has been down-evolved at LO 
using again the HOPPET code ^571] . The solid line is the result for the SSA evolved at LO to the experimental scale 
Q2 = 2.5 GeV2. 

The agreement with data is satisfactory but for the unfavoured channel of n~ emission, probably because of the 
approximation introduced into the description of the Collins function for this case [5^. The effect of DGLAP evolution 
is not large, suggesting that there are compensations between the numerator and denominator of the SSA. 

In Fig. [HI we show the x (above) and z (below) dependences of the weighted SSA sin{<ph-<l>s) ^^.^^ (|T^ for 
both TT"*" (left) and tt~ emission (right) in the same SIDIS kinematics of HERMES as in the previous figures (Sivers 
effect [i^l). Experimental data are taken from Ref. [l^. The dashed line represents the result of the weighted SSA 
when calculating it at the model scale Qq = 0.3 GeV^ by down-evolving the Di of Ref. [55] at LO using again the 
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FIG. 9: The weighted single-spin asymmetry A^j; -^s) ^q. ^ in the SIDIS kinematics of HERMES (Sivers effect HI]) 
with emission of tt"*" (left) and 7r~ (right), as a function of x (above) and z (below). Experimental data from Ref. [l^. Dashed 
and solid lines with the same notations as in previous figure. 



HOPPET code [531 . The solid line is the result for the weighted SSA evolved at LO to the experimental scale = 2.5 
GeV2. 

The agreement between our results and the data is good for the x dependence, while a discrepancy is evident for the 
z distribution in the tt"*" channel. As already anticipated at the beginning of Sec. lIVl the lacking of a strong sea-quark 
contribution at small x (in our model, it is generated only by radiative corrections) may deplete the denominator 
of the SSA once integrated upon x, and produce the observed enhancement upon data in the z distribution (in the 
numerator, the Sivers function is less affected by sea quarks). 

In Fig. \W\ we show our predictions for the x dependence of the weighted SSA A';^^ sin{4>H+4>s) ^ (Collins 

effect), A^^ ^ ^g.^g^g ggg^^)^ ^Qt cosiH-4's) reaction n^(e, eV^). The 

measurement was recently performed at Hall A of JLab by the E06-010 and E06-011 collaborations [24|, using a 6 
GeV energy beam and a transversely polarized '^He (effective neutron) target in SIDIS kinematics with cm energy 
squared s = 12.14 GeV^, average scale (Q^) = 2.2 GeV^ and the following experimental cuts: 



0.13 < x < 0.4 , 0.46 < z< 0.59 , 

1.3 5.4 



Vmin 

(x) < y < 0.86 , yminix) = max 



0.68, 



x{s ~ NP) ' {l-x){s- M2) 



(97) 



In each panel of Fig. [101 the solid (dashed) line represents the result for the emission of 7r+ (tt^) in the above 
kinematics. All components of the weighted SSA have been evolved at LO to the experimental average scale = 2.2 
GeV^ using the HOPPET code [57||. The LO evolution kernels for the transversity and the first Pt moment of the 

Sivers function are the same as those ones used in Fig. [SI and [HJ respectively. The evolution of pj^^ is assumed at LO 
to be the same as the one of the helicity distribution gi . 



C. The SSA in Drell-Yan 



As there are no data for weighted SSA in Drell-Yan collisions, we will show our predictions for cases of interest in 
view of future experiments. At FAIR (GSI), the PAX collaboration is planning to measure the fully polarized Drell- 
Yan process with antiprotons in order to perform a self-consistent extraction of the transversity distribution [25.-29I] . 
Moreover, by simply switching on and off the transverse polarization of only one hadron, from the combination of the 
various cross sections it is possible to extract hi also together with the Boer-Mulders function [see Eqs. (1511 [5^ ]. 
Finally, using the same process the COMPASS collaboration is planning to extract the Sivers function using a high 
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FIG. 10: The weighted single-spin asymmetries A'^'^ ^in(<l>h+'t>s) ^ (CoUins effect), A'^'^ <>i'^(<l>h-i>s) ^ (Sivers 

effect), co^i'Ph-t's) of ([14]), as functions of a; in the SIDIS kinematics at JLab with a 6 GeV beam energy for the 

emission of n'^ (sohd hne) and tt^ (dashed line) from a transversely polarized neutron [2^ . 



energetic pion beam on a transversely polarized proton target at CERN [30|-l33j , in order to directly test its predicted 
non-universal behaviour 1341 1 . 



Arr/aTT 
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FIG. 11: The double-spin asymmetry Att/o-tt of Eq. (|25|) as a function of xf for the p^p^ fi'^ fi~ X process at the PAX 
kinematics (see text). Solid, dashed, dot-dashed, and dotted lines correspond to cm. energies squared s = 30, 45, 80, 200 
GeV^, respectively. 



In Fig. 1111 we show our prediction for the x p dependence of the double-spin asymmetry Att /o,tt from Eq. (|25p 
for the p^p^ — >■ fi^ iJL~X process that could be studied at PAX with a fully transversely polarized antiproton beam 
(|S'it| = 1) colliding on a fully (|S^2r| — 1) transversely polarized proton target [2^[23|. Conventions are according 
to Eq. (P^ . Dimuon invariant masses are summed in the range 2 < M < 3 GeV, below the J /'ip resonance. Solid, 
dashed, dot-dashed, and dotted lines correspond to cm. energies squared s = 30, 45, 80, 200 GeV^, respectively. All 
the parton distributions entering the asymmetry are evolved at LO to each = inside the integration range. 

Despite the fact that Att is roughly the "squared" of the transversity distribution, it does change sign for some 
xp depending on the value of s. As it is evident from Eq. (j22p . for different s the same xp and M probe different 
xi and X2, and, in particular, those where one of the two model transversities for quark up changes sign (see Fig. 8 
of Ref. [iO])- Consequently, the asymmetry of Eq. also changes sign. This feature contrasts with other results 
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in the literature [13, [H, [tI, [t^. Moreover, we also see a monotonically decreasing trend in xp (and also in y), 
but our asymmetry is as large as 10% in modulus at most. This is probably due to cancellations that occur when 
summing upon the invariant mass M, or equivalently when integrating upon specific portions of the {xi,X2) phase 
space, where the product hi{xi) /i"(a;2) in Att repeatedly changes sign [and is comparable with, or bigger than, the 
positive hf{xi) hf{x2)]- For both reasons, a measure of Att is highly desirable because it would also clarify the issue. 
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FIG. 12: The weighted single-spin asymmetry A^Jj, '''^'^'^^^^ /utt of Eq. (|32|) for the pp^ — >■ /i'^ jj, X process as a function of 
sf, in the same kinematics, conventions and notations, as in the previous figure. 

In Fig. ll2[ we show our prediction for the xp dependence of the weig hted SSA i^V'"^^^'^^' VflTT from Eq. m for 
the pp^ — iX^ tx~ X process that could be studied at PAX with antiproton beams colliding on a transversely polarized 
proton target [2^ |43|. Conventions and notations are the same as in the previous figure. Evolution effects are 
included for all partonic densities; as already explained at the end of Sec. Ill Bl the LO evolution equations of the 
chiral-odd /i^^^^ are assumed to be the same as those of the chiral-odd transversity hi. The outcome in Fig. [T^ already 
takes into account the predicted sign change of T-odd TMDs when going from SIDIS to Drell-Yan collisions [s^l • 

Again, for different s the same xp and M probe different xi and X2, in particular the range where the hi{x2) 
changes sign. Consequently, the asymmetry of Eq. (j32l) also changes sign, the turning point depending on s. For the 
very same reason, cancellations probably occur when summing upon the invariant mass M, which suppress the size 

of the asymmetry at positive xp. On the contrary, the size of A'^j,^^'^^'^^'^^^^ /gtt is significant at negative xp. . 
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FIG. 13: The weighted single-spin asymmetry A^Jj, ^'^ of Eq. (ISlf) as a function of a;^ for the vrp^ — >■ /i"*"/! X process at 
COMPASS (see text). Upper solid (lower dashed) line for vr" (tt"''). 

In Fig. 1131 we show our prediction for the x p dependence of the weig hted SSA f^ojn Eq. (EH) for the 

7rp1~ — 7> ^'^ jjL" X process that could be studied at COMPASS with a 160 GeV pion beam colliding on a transversely 
polarized proton target [s^, HH [s^ . Dimuon invariant masses are summed in the safe range 4 < M < 9 GeV between 
the J/V' and the T resonances; since s = 300 GeV^, the invariant r w 0.16 is in the valence region. Other kinematic 
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onventions are the same as in the previous figure. Again, the result already takes into account the predicted sign 
change of T-odd TMDs when going from SIDIS to Drell-Yan collisions ^33| . 

Following Ref. [131, we approximate the A'^rp^™'"'^ of Eq. ([31]) by considering only the dominant valence con- 
tribution to the TT — p collision. Therefore, for 7r~ [t^'^) we retain only the u — u [d — d) annihilation and the final 
weighted SSA does no longer depend on the partonic densities in the pion. Evolution effects are included at LO along 
the lines described at the end of Sec. Ill A] The upper solid (lower dashed) line is our result for tt^ (tt^). 

The size of the asymmetry should allow for an unambiguous test of the prediction about the above mentioned sign 
change of the Sivers function with respect to its extraction from SIDIS [31l] . 

V. CONCLUSIONS 

In this paper, we have analytically calculated azimuthal (spin) asymmetries in lepton-nucleon semi-inclusive deep- 
inelastic scattering (SIDIS) and in Drell-Yan hadronic collisions, using the leading-twist transverse-momentum de- 
pendent distributions (TMDs) obtained in the diquark spectator model of the nucleon from Ref. [lo| . 

We have included evolution effects at leading order (LO) in as by implementing in the HOPPET code of Ref. [s^l 
the kernel for chiral-odd objects. While DGLAP equations are well known for the momentum fi{x), helicity gi{x)^ 
and transvcrsity hi{x) distributions, we have included evolution effects in an approximate way also for the first 

moments of the Sivers f^r^^\x), g|^(a:), and Boer-Mulders hi'^^\x) functions, as described in Sec. IIVI Moreover, 
since in our model the first moment of naive T-odd densities depends linearly upon Us at the model scale, a 
consistent treatment of evolution effects requires to determine itself from the renormalization group equations. 
Therefore, with respect to Ref. [l3| we have modified its value at the hadronic model scale QqI instead of using an 
ad-hoc nonperturbative input as a free parameter, we have computed as((3o) from the LO renormalization group 
equations. The net effect is a resizing of all T-odd functions by a global constant. The comparison with available 
parametrizations [59l . IbsI l66l | is encouraging for the Sivers function, but requires more investigations in the case of the 
Boer-Mulders one. 

We have tested our model TMDs by considering weighted single-spin asymmetries (SSA) in both SIDIS and Drell- 
Yan collisions. Data for weighted SSA are scarce, with low statistics, and still preliminary [H, H^]. But, from the 
theoretical side they are preferable than unweighted SSA because in a model independent way their final expressions 
get automatically factorized involving simple products of "collinear" objects [16| - parton distribution functions 
(PDFs) or Pt moments of TMDs (and analogously for fragmentation functions). As such, weighted SSA can be 
always calculated analytically in our model. Incidentally, the choice of considering evolution effects only at LO is 
dictated by the consistency with the expressions for the weighted SSA, which are known at LO. 

As for SIDIS, we have considered the asymmetry All with longitudinally polarized protons [lH, and the Collins 
and the Sivers effects [T5||, all measured by the HERMES collaboration; we have made predictions for the Collins 
and the Sivers effects, as well as for A^t, in the kinematic conditions recently explored in the E06-010 and E06-011 
experiments in Hall A at JLab [IJl, using the 6 GeV energy beam and the transversely polarized ^He effective neutron 
target. For the unpolarized fragmentation function Di, we have adopted the parametrization of Ref. ^5^, while we 
have used the analytic expression of the Collins function from Ref. [56.] , obtained using a spectator approach similar 
to the present framework. Overall, the comparison with experimental data is satisfactory. For the Collins effect, in 
some cases there are discrepancies that probably can be traced back to the assumptions made in the description of 
the unfavoured channels of the Collins function [56[ . 

Since there are no data for weighted SSA in Drell-Yan collisions, we have made predictions for cases of interest in 
view of future experiments. At FAIR (GSI), the PAX collaboration is planning to measure the fully polarized Drell- 
Yan process with antiprotons in order to perform a self-consistent extraction of the transversity distribution [25-29]; 
we have presented our predictions for Att in various kinematic configurations. Interestingly, despite the fact that 
this asymmetry is approximately the "squared" of the transversity distribution hi, our Att does change sign in some 
portions of the phase space. Briefiy, the kinematics of the parton-antiparton annihilation explores different ran ges in 
the X dependence of the two involved hi, which can separately become negative in our model (see Fig. 8 of Ref. [lO|'). 
This feature contrasts with other results in the literature [13, [H, [tI, [tII ; it would be highly desirable to measure Att 
in order to clarify this issue. 

By simultaneously considering Drell-Yan collisions where both hadrons are unpolarized or only one is transversely 
polarized, it is possible to measure a weighted SSA that leads to the extraction of hi in combination with the Boer- 
Mulders function hi - We have presented our predictions again in the kinematic regime that the PAX collaboration 
could explore at FAIR (GSI). The same previous comment on the sign of hi applies here too. 

Finally, using the same process (single-polarized Drell-Yan collision) the COMPASS collaboration is planning to 
extract the Sivers function using a high energetic pion beam on a transversely polarized proton target at CERN |30l - 
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[33|, in order to directly test its predicted non- universal behaviour (34j. We have shown our predictions for both tt 
and TT"*" collisions, assuming that the elementary annihilation is driven by the dominant valence contributions. 



Acknowledgments 

This work was partially supported by the Research Infrastructure Integrated Activity "Study of Strongly Interacting 
Matter" (acronym IIadronPhysics2, grant agreement n. 227431) under the 7*'' Framework Programme of the Euro- 
pean Community, and by the Italian MIUR through the PRIN 2008EKLACK "Structure of the nucleon: transverse 
momentum, transverse spin and orbital angular momentum" . 



[1] J. C. Collins and D. E. Soper, Nucl. Phys. B194, 445 (1982). 

[2] S. J. Brodsky, D. S. Hwang, and I. Schmidt, Phys. Lett. B530, 99 (2002), hep-ph/0201296. 

[3] S. J. Brodsky, D. S. Hwang, and I. Schmidt, Nucl. Phys. B642, 344 (2002), hep-ph/0206259. 

[4] M. Burkardt and D. S. Hwang, Phys. Rev. D69, 074032 (2004), hep-ph/0309072. 

[5] K. Goeke, S. Meissner, A. Metz, and M. Schlegel, Phys. Lett. B637, 241 (2006), hep-ph/0601133. 

[6] Z. Lu and I. Schmidt, Phys. Rev. D75, 073008 (2007), hep-ph/0611158. 

[7] S. Meissner, A. Metz, and K. Goeke, Phys. Rev. D76, 034002 (2007), hep-ph/0703176. 

[8] J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Rev. D76, 074029 (2007), arXiv:0706.1196 [hep-ph]. 

[9] M. Burkardt and B. Hannafious, Phys. Lett. B658, 130 (2008), arXiv:0705.1573 [hep-ph]. 
[10] A. Bacchetta, F. Conti, and M. Radici, Phys. Rev. D78, 074010 (2008), 0807.0323. 
[11] S. Chekanov at al. (ZEUS), Phys. Rev. D67, 012007 (2003), hep-ex/0208023. 

[12] M. Gliick, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D63, 094005 (2001), hep-ph/0011215. 
[13] A. Airapetian et al. (HERMES), Phys. Rev. Lett. 94, 012002 (2005), hep-ex/0408013. 
[14] A. Bressan (COMPASS) (2009), 0907.5508. 

[15] R. Seidl (HERMES) (2004), prepared for 12th International Workshop on Deep Inelastic Scattering (DIS 2004), Strbske 

Pleso, Slovakia, 14-18 Apr 2004. 
[16] P. J. Mulders and R. D. Tangerman, Nucl. Phys. B461, 197 (1996), erratum-ibid. B484 (1997) 538, hep-ph/9510301. 
[17] M. Guagnelli, A. Bacchetta, and M. Radici (in preparation). 

[18] S. Boffi, A. V. Efremov, B. Pasquini, and P. Schweitzer, Phys. Rev. D79, 094012 (2009), 0903.1271. 

[19] Z.-B. Kang and J.-W. Qiu, Phys. Rev. D79, 016003 (2009), 0811.3101. 

[20] W. Vogelsang and F. Yuan, Phys. Rev. D79, 094010 (2009), 0904.0410. 

[21] V. M. Braun, A. N. Manashov, and B. Pirnay Phys. Rev. D80, 114002 (2009), 0909.3410. 

[22] A. Airapetian et al. (HERMES), Phys. Rev. D71, 012003 (2005), hep-ex/0407032. 

[23] M. Alekseev et al. (COMPASS), Phys. Lett. B680, 217 (2009), 0905.2828. 

[24] JLab experiment E06-010/E06-011, J. -P. Chen, E. Cisbani, H. Gao, X. Jiang, J.-C. Peng, spokespersons. 
[25] V. Barone et al. (PAX) (2005), hep-ex/0505054. 

[26] M. Maggiora et al. (ASSIA), Czech. J. Phys. 55, A75 (2005), proceedings of the ASI Conference on Symmetries and Spin 

(Spin-Praha 2004), Praha, 5-10 July 2004., hep-ex/0504011. 
[27] A. V. Efremov, K. Goeke, and P. Schweitzer, Eur. Phys. J. C35, 207 (2004), hep-ph/0403124. 
[28] M. Anselmino, V. Barone, A. Drago, and N. N. Nikolaev, Phys. Lett. B594, 97 (2004), hep-ph/0403114. 
[29] A. Bianconi and M. Radici, Phys. Rev. D72, 074013 (2005), hep-ph/0504261. 

[30] A. V. Efremov, K. Goeke, S. Menzel, A. Metz, and P. Schweitzer, Phys. Lett. B612, 233 (2005), hep-ph/0412353. 
[31] A. Bianconi and M. Radici, Phys. Rev. D73, 114002 (2006), hep-ph/0602103. 
[32] A. Bianconi and M. Radici, Phys. Rev. D73, 034018 (2006), hep-ph/0512091. 
[33] M. Anselmino et al., Phys. Rev. D79, 054010 (2009), 0901.3078. 
[34] J. C. Collins, Phys. Lett. B536, 43 (2002), hep-ph/0204004. 

[35] A. Bacchetta, M. Diehl, K. Goeke, A. Metz, P. J. Mulders, and M. Schlegel, JHEP 02, 093 (2007), hep-ph/0611265. 
[36] A. Bacchetta, D. Boer, M. Diehl, and P. J. Mulders, JHEP 08, 023 (2008), 0803.0227. 
[37] D. Boer and P. J. Mulders, Phys. Rev. D57, 5780 (1998), hep-ph/9711485. 
[38] J. Zhou, F. Yuan, and Z.-T. Liang (2009), 0909.2238. 

[39] V. Barone, A. Drago, and P. G. Ratcliffe, Phys. Rept. 359, 1 (2002), hep-ph/0104283. 
[40] M. Stratmann and W. Vogelsang, Phys. Rev. D65, 057502 (2002), hep-ph/0108241. 
[41] J. C. Collins, D. E. Soper, and G. Sterman, Nucl. Phys. B250, 199 (1985). 
[42] J. C. CoUins and D. E. Soper, Phys. Rev. D16, 2219 (1977). 
[43] S. Arnold, A. Metz, and M. Schlegel (2008), 0809.2262. 

[44] R. D. Tangerman and P. J. Mulders, Phys. Rev. D51, 3357 (1995), hep-ph/9403227. 

[45] D. Boer, Phys. Rev. D60, 014012 (1999), hep-ph/9902255. 

[46] A. Bianconi and M. Radici, Phys. Rev. D71, 074014 (2005), hep-ph/0412368. 



21 



[47] M. Radici, F. Conti, A. Bacchetta, and A. Bianconi (2007), arXiv:0708.0232 [hep-ph], URL 
Ihttp : / / www . slac . stanf ord . edu/ spires/f ind/hep/ www?eprint=arXiv : 0708 . 0232 
R. Jakob, R J. Mulders, and J. Rodrigues, Nucl. Phys. A626, 937 (1997), liep-ph/9704335. 
H. Avakian, A. V. Efremov, P. Schweitzer, and F. Yuan, Phys. Rev. D78, 114024 (2008), 0805.3355. 
A. V. Efremov, P. Schweitzer, O. V. Teryaev, and P. Zavada (2008), 0812.3246. 
J. She, J. Zhu, and B.-Q. Ma, Phys. Rev. D79, 054008 (2009), 0902.3718. 
D. W. Sivers, Phys. Rev. D41, 83 (1990). 

D. Boer, P. J. Mulders, and F. Pijlman, Nucl. Phys. B667, 201 (2003), hep-ph/0303034. 

A. Bacchetta, M. Boglione, A. Henneman, and P. J. Mulders, Phys. Rev. Lett. 85, 712 (2000), hep-ph/9912490. 
D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. D75, 114010 (2007), hep-ph/0703242. 
A. Bacchetta, L. P. Gamberg, G. R. Goldstein, and A. Mukherjee, Phys. Lett. B659, 234 (2008), 0707.3372. 
G. P. Salam and J. Rojo, Comput. Phys. Commun. 180, 120 (2009), 0804.3755. 
M. Gliick, E. Reya, and A. Vogt, Eur. Phys. J. C5, 461 (1998), hep-ph/9806404. 
M. Anselmino et al., Eur. Phys. J. A39, 89 (2009), 0805.2677. 
J. C. Collins et al. (2005), hep-ph/0510342. 

M. Diefenthaler (HERMES) (2007), proceedings of the 15th International Workshop on Deep Inelastic Scattering (DIS 
2007), Munich, Germany, 16 - 20 Apr 2007., 0706.2242. 
A. Martin (COMPASS), Czech. J. Phys. 56, F33 (2006), hep-ex/0702002. 
M. Alekseev et al. (COMPASS), Phys. Lett. B673, 127 (2009), 0802.2160. 
M. Gockeler et al. (QCDSF), Phys. Rev. Lett. 98, 222001 (2007), hep-lat/0612032. 
V. Barone, S. Melis, and A. Prokudin (2009), 0912.5194. 
Z. Lu and I. Schmidt (2009), 0912.2031. 
W. Kafer (COMPASS) (2008), 0808.0114. 

A. Bressan (COMPASS) (2009), 0907.5511. 

F. Giordano and R. Lamb (On behalf of the HERMES), AIP Conf. Proc. 1149, 423 (2009), 0901.2438. 
L. Y. Zhu et al. (FNAL E866/NuSea), Phys. Rev. Lett. 102, 182001 (2009), 0811.4589. 
L. P. Gamberg, G. R. Goldstein, and M. Schlegel, Phys. Rev. D77, 094016 (2008), 0708.0324. 
F. Yuan, Phys. Lett. B575, 45 (2003), hep-ph/0308157. 
J. C. Collins, Nucl. Phys. B396, 161 (1993), hep-ph/9208213. 

V. Barone, A. Cafarella, C. Coriano', M. Guzzi, and P. Ratcliffe, Phys. Lett. B639, 483 (2006), hep-ph/0512121. 

B. Pasquini, M. Pincetti, and S. Both, Phys. Rev. D76, 034020 (2007), hep-ph/0612094. 



